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Abstract. The harmonic Fibonacci chain, which is one of a quasiperiodic 
chain constructed with a recursion relation, has a singular continuous frequency- 
spectrum and critical eigenstates. The validity of the Fourier law is examined 
for the harmonic Fibonacci chain with stochastic heat baths at both ends by 
investigating the system size N dependence of the heat current J and the local 
temperature distribution. It is shown that J depends on Af as J ~ (InA^)"^ 
and the local temperature strongly oscillates along the chain. These results 
indicate that the Fourier law does not hold on the harmonic Fibonacci chain. 
Furthermore the local temperature exhibits two different distribution according to 
the generation of the Fibonacci chain, i.e., the local temperature distribution does 
not have a definite form in the thermodynamic limit. The relations between N- 
dependence of J and the frequency-spectrum, and between the local temperature 
and critical eigenstates are discussed. 



PACS numbers: 44.10.-|-i, 61.44.Br, 05.40.-a, 05.60.-k 



Many studies in these decades have shown that arbitrarily defined one-dimensional 
(ID) systems of interacting particles do not exhibit normal thermal transport 
properties, i.e., the Fourier law does not hold on such systems [|l|, H, ||, |[ ||. For 
the steady-state of the homogeneous ID chain of system size N, the Fourier law, 
J = —kS/T with thermal conductivity k, indicates that the heat current depends 
on the system size as J ~ 1/N and that the temperature gradient VT is constant 
along the chain. Rieder et al Q have shown for a ID chain of equal mass particles 
interacting with identical harmonic potential that the heat current is independent of 
the system size. They also have obtained the local temperature distribution. The local 
temperature behaves in unphysical way: the temperature takes a constant value in the 
bulk. Furthermore, near the end of the chain, the temperature decreases as we move in 
the direction of the hotter heat bath, and rise only at the end particle in contact with 
the heat bath; the temperature exhibits corresponding behavior at the other end of 
the chain. Casher and Lebowitz |2] have shown for the same model but with random 
mass distribution that J ^ N^^''^. For the same random mass distribution model 
but with different type of boundary conditions, Rubin and Greer ^ have obtained 
the result J ~ iV^^/^. These system size dependence of J for periodic or disordered 
chains may be attributed to the localization property of eigenstates on these chains. 
Since the eigenstates are extended in periodic chains, the ballistic energy transport 
of extended eigenstates results in the constant heat current of the periodic chains. 
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In contrast, for the disordered chains, the decrease in heat current with increasing 
system size is caused from the locahzed eigenstates, which can not transport energy 
over the length of the system. The extended and locahzed eigenstates correspond to 
continuous and pure-point spectra, respectively. It is shown that, for the Casher and 
Lebowitz type chain and heat bath, the thermal conductivity k diverges as the system 
size increases if the spectrum contains absolutely continuous part 

Some of quasiperiodic systems have a Cantor set-like spectrum, i.e., a singular 
continuous spectrum [|[ ^, ||, |[ |l^ . Then, on such systems, eigenstates show power 
law decay; such neither extended nor exponentially localized states are called critical 
states For example, a harmonic chain with quasiperiodically arranged spring 
constants and/or mass of particles has a singular continuous spectrum and critical 
eigenstates 0, |[ ||] ■ We may thus expect that such quasiperiodic systems show exotic 
heat transport properties compared with periodic or disordered systems. 

In the present paper we investigate the anomaly of heat transport phenomena on 
harmonic Fibonacci chains, which is one of discrete quasiperiodic and self-similar ID 
lattices. We focus the anomaly resulting from the spectral properties of the Fibonacci 
chain. In order to check the validity of the Fourier law on the Fibonacci chain, we 
investigate the system size dependence of the heat current J. Although Macia have 
already studied the thermal conductivity k of the harmonic Fibonacci chain, he has 
not checked the system size dependence of k and thus the validity of the Fourier law is 
not clear yet. Our results show that the heat current behaves as J ~ (InTV)"-'^, which is 
in contrast with that of periodic or disordered chain. We discuss the fact that the total 
bandwidth of the phonon spectrum of the Fibonacci chain has similar dependence 
to J. We also calculate the local temperature distribution on the Fibonacci chain; it 
seems not to converge to a definite form even in the thermodynamic limit. We relate 
{Ti} with the critical eigenstates of the Fibonacci chain. 

The harmonic Fibonacci chain which we consider is a ID chain of N particles; 
each particle interact with its neighbouring particles with equal spring constant k. We 
make the sequence of mass of particles {mi|i = 1, ...,N,mi = rria or m/j} according 
to the Fibonacci sequence. The Fibonacci sequence of the n-th generation L„, which 
consists of two kinds of components rua and mp, is constructed by the recursion 
relation: L„ = L„_iL„_2, Lq = m^, and Li — nia- Then the system size of the n-th 
generation Fibonacci sequence is the Fibonacci number Fn , which obeys the recursion 
relation Fn = Fn-i + Fn-2, — 1, and Fi = 1. The Fibonacci number F„ is 
asymptotically behaves as F„ ~ r" with golden ratio r = {V5 + l)/2. We can obtain 
the asymptotic properties of the Fibonacci chain in the limit of A'^ — > oo by concerning 
the infinite-generation limit n — > oo. We set both the spring constant k and the mass 
unity; we calculate the heat current and the temperature distribution varying with 
the mass rUa- 

We consider the following Langevin equations for particles of the chain with 
stochastic heat baths at both ends: 

mixi = -2x1 + X2- 7^1 + rji^it), 

rriiXi^ -2xi+ Xi^i+ Xi+i, (1) 

rriNXN = -~'2xN + xn~i — JXn + 

where Xi are the displacements of the particles from their equilibrium positions; 7 is 
the friction constant; t^l and tjb, are the random forces caused from left and right heat 
baths, respectively. We choose the random forces the white noises, i.e., the Fourier 
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transforms of the random force obey 

(?7l(w)?7l(w')> = 4:Tr'yT-LS{uj + u'), 
{r]R{uj)ri-R{uj')) = ATrjTjxd{uj +uj'), 

where angular bracket is the average over the random force; Tl and Tr are the 
temperatures of the left and right heat baths, respectively. 
The energy current Ji from the site i to the site i + 1 is 

Ji ^ ^{{ii+i + Xi){xi+i - Xi)). (3) 

In the steady state, the energy current does not depend on the site; then from the 
Langevin equation (|lj), 

J-AtT-^^ (4) 

-^-"^'Jo ^ Idetyp' 

where AT = Tl - Tr and F = $ - uj'^M + iujT with ^ij = + 2Sij - Sij+i, 

Mij = rriiSij, and = 7((5,;,i + 5i^N)5i,j- We may write 

det Y = Di^N - 7^w^£)2,Ar_i + i^uj{D2,N + -Di,jv-i), (5) 

where Di j denotes the determinant of the sub-matrix of (f> — cu'^M which begins with 
i-th row and column and ends with j-th row and column Q . 
We define the local temperature at site i 

where H is the Hamiltonian of the system. We note that the equipartition relation 
Ti = (xidH/dxi) holds for any mass distribution (for periodic case, the equipartition 
relation has proven in [Q). For our model we can write the explicit form of the local 
temperature as 

Z^Li-i + i7LjT>2,i-i ^ 

• UJ 

-oo 

In calculating the heat current, it is convenient to express the determinant Dij 
in terms of a transfer matrix as 

jF'^' Sn'^'^' ) = M,M,^i ■ ■ ■ M,+iM,, (8) 

where the matrix Mi{uj) is the unimodular transfer matrix of site i defined as 

M^{UJ) ^[ I Q )■ (9) 

If we rewrite the total transfer matrices of the n-th generation Fibonacci chains as 
Ain = MnMn^i ■ ■ ■ Ml with N ^ Fn, then A1„ obey a recursion relation 

Mn = Mn-2Mn-i, AIq = ^ ^ "^^^^ ) , and Ml - Ah. (10) 

Since we can obtain all the determinants Dij needed in (^) by calculating Mn, the 
recursion relation ( p^ saves a lot of time in calculating J. 

We calculated the heat current and the local temperature by numerical integration 
of (^ and (|^), respectively. As the system size grows larger, and as the mass ratio 
ma/mp becomes farther from unity, the integrands of (^) and (^ vary more rapid; 
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hence we should use finer intervals on numerical integration. For all the following data 
of J and {Tj} we checked that decreasing the intervals of integration by a factor 5 did 
not make significant changes of the results. 

From the numerical calculation, we found that the total energy current J decreases 
with oscillation as the generation n increases (see figure |l|) . The decrease of J for the 
generation n shows power law behaviour J ~ n~°; the exponent is a ~ 1 for sufficiently 
large rria (e.g., a = 1.03 when rria = 3.0 ). This result indicates that it is appropriate 
to discuss the heat current in terms of the generation n rather than the system size 
N. We note that the system size dependence J ^ n'"" ^ (IniV)"" is the remarkable 
difference from the behaviour J ~ of periodic or disordered systems. 

In order to discuss the contribution of eigenstates of the chain to the total heat 
current we calculate, as done in |pT| , cumulative heat current: 



The cumulative heat current Jc(w) is a monotonically increasing function; the 
increment of Jc{uj) in a certain w-region is the contribution of the eigenstates of the 
w-region to the total heat current. For the Fibonacci chain, as shown in figure]^, Jc{uj) 
consists of many slopes and plateaus. Each slope consists of slopes and plateaus in 
nested fashion, and the nested structure becomes finer as the generation grows; i.e., 
Jc(tt') seems to become the devil's staircase in the infinite-generation limit. 

The devil's staircase-like structure of Jc{uj) is very similar to the integrated 
density of states (IDOS) of the Fibonacci chain (see figure ^). The plateaus of 
Jc(w) approximately coincide with the plateaus of the IDOS; i.e., the frequencies 
corresponding to the eigenstates make major contribution to the heat current. 

Figure ^ shows that dominant contribution to J comes from the low frequency 
eigenstates, especially in high generation. We may attribute the high conductivity in 
the low-frequency region to the extended eigenstates. The high frequency eigenstates 
are critical, i.e., they show a power-law decay, while the low frequency eigenstates 
are extended like eigenstates in periodic chains. Correspondingly, the low frequency 
eigenstates may transport the energy ballistically. The ballistic transport can be seen 
from the similarity between Jc{^) of the Fibonacci chain and Jc{^) of a periodic chain 
in low frequency region (see figure |^) . We should note that the low frequency region 
in which the eigenstates are extended become narrower in the limit of infinite system 
size. This is because all the eigenstates become critical in the limit of infinite system 
size, since the continuous part of the spectrum of the Fibonacci chain narrows and 
becomes singular continuous even in low frequency region Q. As a consequence, we 
can conclude that the heat current, unlike that of the periodic chain, vanishes when 
N ^ oo. 

The similarity between Jd^^) and IDOS indicates that the generation dependence 
of J is explained by the spectral properties of the Fibonacci chain. Now, we consider 
the n-th generation rational approximants, which is a periodic chain whose unit cell 
is the n-th generation Fibonacci chain. We denote by cr„ the w^-spectrum of the n-th 
generation rational approximant; (T„ is given as 

{cj2; I tr7W„| < 2.}. The spect rum a„ 
consists of Fn bands; as n increases, the bands are fragmented into pieces and total 
bandwidth W„ decreases. We note that Wn obeys a power law Wn ~ n~° and a' ~ 1 
0. Since the frequency region corresponding to the bandgap does not contribute to 
the heat current, we conclude that the decay of J with increasing n is similar to the 
decay of W„; thus J obeys the power law as Wn does, as shown in figure |l|. 
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Before showing the local temperature {Tj}, we define Fibonacci sites and 
Fibonacci blocks, which seem important to discuss the behavior of {T^}. On the 
n-th generation Fibonacci chain, which consists of sites, the Fibonacci site is the 
Fj-th site with j = 1, 2, .... n — 1 from left or right ends; i.e., the Fibonacci sites are 
numbered Fj or F„ — Fj . We call a region between neighbouring Fibonacci sites the 
Fibonacci block. Since Fj ~ , Fibonacci sites are sparsely distributed around the 
centre of the chain and densely near the ends. If the chain length is normalized to unity, 
the Fibonacci sites are approximately at and 1 — . Thus, as the generation 
increases, the Fibonacci sites distribution becomes denser near the ends but invariant 
around the centre of the chain. Correspondingly, the size of the Fibonacci blocks 
around the centre of the chain stays constant in increasing generation. 

Now we plot {Ti} with Fibonacci sites in figure |^. We can see that the local 
temperature oscillates strongly along the chain; especially {Ti} has significant leaps 
at the Fibonacci sites. Furthermore, the interior of each Fibonacci blocks have a 
bilaterally symmetric temperature distribution. The symmetry is a nontrivial result 
since the whole system with heat baths at both ends is not bilaterally symmetric. We 
note that some of the Fibonacci blocks have a symmetric mass sequence, but not all 
the blocks do. 

Another major characteristic of {Ti} is an oscillatory behavior for increasing 
generation. We plot in figure ^ the local temperature on the Fibonacci chains of 
generation n — 12, 13, 14 and 15. There is an obvious similarity between {Ti} of 
n = 12 and n — 14, and between {Ti} of n = 13 and n — 15; i.e., {Ti} oscillates 
periodically in generation with period two. We checked that the similarity in {Ti} of 
odd- and even-generation chains is maintained throughout the generations from n = 6 
to n = 15, which correspond to the system size from iV = 13 to iV = 987. We may 
check the oscillation quantitatively by comparing the average local temperature of, 
e.g., the centre Fibonacci block of each chain. The average temperature of the centre 
Fibonacci block is distributed between 1.558 and 1.589 for odd-generation chains (from 
n = 7 to n = 15); between 1.645 and 1.706 for even-generation chains (from n = 6 to 
n — 14). The periodic oscillation implies that the local temperature distribution {T,} 
does not converge to a definite distribution even in the limit of infinite-generation, 
i.e., in the thermodynamic limit. 

It is interesting to note that the mass-weighted amplitude of critical eigenstates 
also have leaps at the Fibonacci sites as shown in figure Such behaviour is similar 
to that of the {Ti}; thus we may expect that the characteristic of the {Ti} partly 
owes to the critical eigenstates. Since we do not observe a critical eigenstate with the 
symmetry in the Fibonacci block, we may not attribute the cause of the symmetry of 
{Ti} in Fibonacci sites to the properties of a single critical eigenstate. 

In summary, we have shown that the Fourier law is invalid for the harmonic 
Fibonacci chain as for periodic or disordered harmonic chains. The heat current 
J depends on the system size N as J ^ (IniV)"-^; the anomalous dependence is 
remarkably different from the periodic or disordered case. The N dependence of J is 
very similar to the N dependence of the total bandwidth of the frequency spectrum. 
Such similarity is the consequence of the fact that only the lattice vibration with 
frequency within the frequency-bands can transport the heat over the length of the 
chain. We have also obtained the local temperature. The local temperature {Ti} 
strongly oscillates and does not show monotonic change. Both {Ti} and the mass- 
weighted amplitude of critical eigenstates have significant leaps at Fibonacci sites; 
this fact indicates the close relation between {Ti} and eigenstates. We have concluded 
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that {Ti} does not converge to a definite form in the hmit of infinite system size since 
{Ti} exhibits two different form depending on the parity of the generation. 
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Figure 1. Log- log plot of the generation n versus the heat current J (circle) 
and the total bandwidth Wn (diamond). Lines are included to guide the eye. The 
masses of particles in the Fibonacci chain are met = 1.9 and = 1.0. The total 
bandwidth Wn is normalized by the maximum eigenfrequency. 
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Figure 2. Full lines indicate the cumulative heat currents Jc{io). The lower two 
full lines are Jcii^) of 10th (upper line) and 20th (lower line) generation Fibonacci 
chains with nia = 1-9 and = 1.0. The upper most smooth full line is Jc{oj) 
of a periodic chain with identical particles; the mass of particles is the average of 
particle in the above two Fibonacci chains. Dashed line is the IDOS of the 15th 
generation Fibonacci chain. Although there are eigenstates in the high-frequency 
region (up to a; ~ 1.7), we plot only low- and middle-frequency region because 
high frequency states contribute little to Jc(w). 
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Figure 3. Upper full line is local temperature {Ti} on the 15th generation 
Fibonacci chain with ma = 1-9 and = 1.0. Triangles indicate the Fibonacci 
sites: upward and downward triangles correspond to Fj-th and (F15 — Fj)-th sites, 
respectively. Lower full line is miX^{uj) of a critical eigenstate corresponding to 
the maximum eigenfrequency. 
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Figure 4. Local temperature on Fibonacci chain of (a) generation n = 12 
(number of particles N = 233), (b) n = 13 (TV = 377), (c) n = 14 {N = 610), and 
(d) n = 15 {N = 987). The Fibonacci chains consist of particles of rria = 1-9 and 
mg = 1.0. Each chain length is normalized to unity. 



